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We present the first principle theory of the existence and stability of TE and TM spatial solitons in 
a subwavelength periodic semiconductor-dielectric structure. We have found that for the wavelength 
1550nm and the interface separation close to and less than lOOnm the band structure of the linear 
TE and TM waves becomes similar to the band structure of a homogeneous medium. The properties 
of TE solitons change accordingly. However, the transverse profiles of the TM solitons continue to 
reflect the subwavelength geometry of the structure and develop dominant intensity peaks inside the 
low index dielectric slots. Our stability analysis based on the linearized Maxwell equations indicates 
that the nonlinear TM waves can be approximated as the evanescently coupled modes of the slot 
waveguides with the low index dielectric core and the high index semiconductor cladding. Transition 
to the regime where the slot waveguides start to determine properties of TM waves is associated 
with the so called Brewster condition. 

PACS numbers: 42.65.Tg,78.67.Pt,42.70.Nq 



I. INTRODUCTION 

Recent progress in the fabrication of nanostructures 
for photonics applications has stimulated research into 
light trapping and guiding on the subwavelength scale 
[H [2] [3. Surface plasmon polaritons tightly confined 
to the metal-dielectric interfaces have been at the focus 
of recent efforts in this direction [H |5l |6l [7]. The re- 
search into plasmons has included studies of their inter- 
action with periodically structured metals, see, e.g. [8 , 
which have been recently extended to soliton structures 
in periodic arrays of metalic slot waveguides filled with 
a nonlinear dielectric [9 . The latter work has been a 
significant and conceptual departure from research into 
optical solitons in nonlinear waveguides arrays coupled 
by evanescent waves [1^. The metal dielectric interfaces 
in [9 are separated by distances much smaller than the 
wavelength, so that the concept of the coupling induced 
by the evanescent waves [lOl [11] becomes largely irrele- 
vant. New effects can be expected in this regime, which 
are still waiting to be explored. 

Semiconductors also can be used for subwavelength 
guiding. In particular, silicon photonic wires have been 
recently promoted as promising and close to practical ap- 
plications building blocks of photonic chips, where non- 
linear and soliton effects have been already extensively 
researched [l2]. The large refractive index of silicon 
(n 3.5) allows for tight light confinement by the con- 
ventional total internal reflection mechanism, giving the 
simultaneous advantages of strong ultrafast Kerr non- 
linearity (n2 4 X 10~^^m^/W), controlled dispersion 
and manageable losses. Losses are a particular problem 
for plasmons, suggesting that their nonlinear functional- 
ity is likely to become more viable if gain is introduced 
into the dielectric [13l [14]. Two photon and free career 
induced absorptions are traditionally thought of as ham- 
pering the attractiveness of silicon for nonlinear applica- 
tions, but these often can be lived with [15 , 16 or man- 
aged, e.g. by the electrically removing free carriers from 
the waveguide core |T2j [TT] [18]. There are also other 



highly nonlinear semiconductors and glasses which can 
be useful for various on-chip applications. In particular, 
many of the soliton experiments in planar waveguide ar- 
rays have been performed using doped GaAs (n = 3.47, 
n2 = 3.3 X 10-^^mVW) structures [E]. 

While an isolated silicon or GaAs photonic wire con- 
fines light within an area of the order of the wavelength 
inside the material squared (A^^^/n^), bringing two wires 
together with a separation of few tens of nanometers pro- 
duces a strong intensity peak in-between the wires^ with 
the field predominantly polarised perpendicular to the 
interface between the wires (TM-modes). This type of 
waveguides is called slot waveguides and they provide an 
elegant method of focusing light into a subwavelength 
region [20 . One can ask a question about the exis- 
tence of solitary waves having subwavelength dimensions 
in arrays of the semiconductor-dielectric slot waveguides. 
The losses in such arrays are expected to be few dB/cm 
[2T| \22\ [23] . Measurements of the transverse profile of the 
slot mode have been reported in [24 , while applications 
of the silicon slot waveguides have been demonstrated for 
frequency conversion and wavelength division multiplex- 
ing as well as for the design of high Q resonators 
[23] and sensors based upon them [26 . 

In this work, we consider an infinite array of semi- 
conductor photonic wires with the wire widths and sep- 
arations taken well below the wavelength. We assume 
that the photon energy is below the bandgap and the 
semiconductor acts essentially as a high index nonlin- 
ear dielectric. We present analysis of linear and soliton 
solutions for both TE and TM polarizations using first 
principle Maxwell equations. The TE solitons undergo a 
smooth transformation from the case, which can be well 
understood in terms of the evanescently coupled waveg- 
uides into the solitons of a quasi-homogeneous medium, 
when the widths of all layers become much less than the 
wavelength. Conversely and centrally to our work, the 
TM solitons, under the same change of geometry, evolve 
into structures with the dominant intensity peaks located 
outside the semiconductor, i.e. in the low index material. 
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FIG. 1: (Color online) One dimensional periodic structure 
consisting of layers of high index semiconductor such as silicon 
(width w) and low-index dielectric such as silica glass (width 
s). 

These peaks become the prevailing features as the separa- 
tion of high index semiconductor waveguides is reduced. 
Thus the TM solitons are not transformed into the soli- 
tons of the quasi-homogeneous medium, which does not 
distinguish between polarizations. Our analysis strongly 
suggests that for subwavelength separations, the TM soli- 
tons can be qualitatively considered as discrete solitons 
composed of the coupled modes of the slot waveguides. 
We also propose and apply a technique to study linear 
stability of solitons within the framework of the linearised 
Maxwell equations. 



II. MAXWELL EQUATIONS AND SOLITON 
EQUATIONS 

The array of slot waveguides we consider below is a 
periodic structure of narrow layers of a high index semi- 
conductor material (material 5, e.g., silicon) embedded 
into a low index dielectric material (material ^, e.g., sil- 
ica glass), see Ref. [21 and Fig. IT] The separations 
of the semiconductor layers vary from 500 to 50nm, and 
the width of the layers is between 220 and 95nm. The 
vacuum wavelength is taken to be Xyac = 1.55/im. Our 
analysis is based on the nonlinear Maxwell equations in 
the 2 dimensional geometry: 

^ ^ ^ ^ ^ k ^ 

\/ X H = -ikceoD, VxE = i — H, (1) 

ceo 

where k = 27r / Xyac ^ c is the speed of light in vacuum, eo is 
the vacuum permittivity, and for electric S and magnetic 
H fields it is assumed that E^H = H-ex.p{—ikct)-\-c.c. 
Light propagation is assumed to be along the z-direction 
and the x-coordinate is perpendicular to the layers, see 
Fig.[l] 

For the TM-polarized modes (such that only Hy^E^ 
and Ez are non-zero) the Maxwell equations are reduced 
to 

d,zE,-d,,E, = -k^D,, (2) 
dz,E,-d:,:,E, = eD,, (3) 
dzHy = ikceoD^. (4) 

While for the TE-polarized modes, only Ey^ and Hz 



are non-zero and the resulting equations are 

dzzEy = -k^Dy, (5) 
k k 

— Hj, = idzEy, — Hz = -ida^Ey. (6) 
ceo ceo 

We only need to solve these equations for the electric 
field components in both the TE and TM cases. The 
constitutive relation is taken as for isotropic materials 

3 = eE+^-X3[\E\^E+^{E ■£)£*], (7) 

where D is the displacement in SI units normalized to 
eo- The above expression for D is an approximation for 
anisotropic semiconductors, but its use is sufficient to 
demonstrate the reality of the effects we are interested in 
and helps to improve the transparency of our results and 
to simplify the complex numerical calculations. We also 
neglect linear and nonlinear absorptions, which influence, 
but do not prevent observation of soliton effects |T5l [16] . 

In our numerical approach, we do not force boundary 
conditions at the interfaces, but instead assume that the 
linear permittivity e and the nonlinear susceptibility xs 
change continuously (but sharply) between their respec- 
tive values for the material s (silicon) and material g 
(silica glass). We model the structure by taking 

e{x) = e^ + ^(e, - eg)Kj{x), (8) 
j 

where Kj{x) = exp{ — [(x — Xj)/wY^} is the array of 
super-gaussian functions, j = 0, ±1, ±2 . . . , = j{s-\-w) 
is the position of the jth semiconductor layer, w is its 
width and s is the side to side separation of the semi- 
conductor layers (see Fig. 1), e = and Ug = 1.44, 
Us = 3.48. xs is linked to the n2 coefficient (measured 
in vr? jW) and found in the tables as X3 = |n2eeoc [27j . 
Thus it is convenient to introduce the function 

7(x) = eocn2(x)e(x), (9) 

where 

n2{x) = n2^g + ^{n2,s - n2^g)Kj{x), (10) 
j 

n2,g = 2.5 • lO-^^mVW, ns,, = 4 • lO-^^m^/W . 

A. Equations for TM solitons 

We seek soliton solutions of Eqs. ([2|, ([3|, ^ in the 
form 

i?. = f{x)e'''''\ = ig{x)e"i''\ (11) 

After some algebra, we find that the real functions / and 
g obey the system of the first order ordinary differential 
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equations 



. _ kdF/dg-q[e' + i{p+g''mf 



<z{e + 7[3/2+5V3]} 
kdF/df 



'/{e + 7[3/2 + 5V3]}^ 



(13) 



where prime indicates first derivative in x and the pa- 
rameter q measures the relative change of the propaga- 
tion constant with respect to its vacuum value k. The 
function F{f^g) is given by 



F = 



2 4 



7 + 



6 



r 

2 



(14) 



F becomes the first integral of Eqs. (12)-(13) in the case 



of a homogeneous medium (when e and 7 are x indepen- 
dent) [28]. 



B. Equation for TE solitons 

Soliton solutions of Eqs. ([5|, ([7| are sought in the form 

Ey = u{x)e'''^' , (15) 

which results in the familiar stationary nonlinear 
Schrodinger equation (NLSE) 



u + r (e - q^)u + k^ju^ = 0. 



(16) 



In a case of the evanescently coupled semiconductor 
waveguides, the above equation is readily transformed 
into a set of coupled mode algebraic equations, giving fa- 
miliar discrete soliton solutions [10 . In the limit when w 
and s are much less then Xyac^ one should expect the soli- 
ton profiles to be close to the ones known from the NLSE 
with constant coefficients. Such qualitative conclusions 
are, however, difficult to make simply by looking at the 
equations for TM solitons. A useful insight into the dif- 
ference between the TM and TE waves can be obtained 
when considering linear modes of the structure. 



III. LINEAR MODES, BAND STRUCTURE 
AND BREWSTER CONDITION 



For 7 = Eqs. (12)-(14) for TM waves are reduced to 
the linear eigenvalue problem 



/■ (17) 



It can be seen that the sharp jumps in e can be compen- 
sated only if the function / itself changes sharply. The 
jump in / is determined by the refractive index contrast, 
and follows from the continuity of at an interface. 



For TE waves the linear equation is 
k^q^u = u" + k'^eu. 



(18) 



Therefore, jumps in e only force the second derivative of 
u to change accordingly, while u itself stays smooth. 

According to the Floquet-Bloch theorem, linear modes 
for TE and TM cases (/ or u) can be represented in the 
form r{x) eyL-^iikskx)^ where r{x^w^s) = r{x) and w-\-s 
is the period. All eigenvalues > are parameterized by 
the Bloch wavenumber /c^, < < 7r/{k{w + s)). In 
the linear case, the problem is tractable analytically with 
exact boundary conditions at the interfaces [29, 30^. The 
discussions in ^3 ^ are focused on the band structure 
and miss important for us features of the linear mode 
profiles, which are highlighted below. 

If we choose w = 220nm and a sufficiently large sep- 
aration between the semiconductor layers s = 500nm, 
then the spectrum ^^(/cb) for both TE and TM modes 
has two bands, see Figs. |2|a) and (b). The ks = and 
= 7r/{k{w + s)) Bloch modes of the top band consist, 
respectively, of the in-phase and out-of-phase modes of 
the individual waveguides, see insets in Figs.[2|a) and (b). 
The TE Bloch modes can be well approximated using the 
well known tight-binding approach [10 , when the Bloch 
mode is represented by a superposition of the modes of 
the individual semiconductor layers. For TM modes, the 
field structure is dominated by jumps at the boundaries. 
However, the field overlap in the dielectric (silica) layers 
still happens via exponentially decaying tails, and so the 
tight binding approximation is justifiable here as well. In 
both cases, the tight binding model will readily reproduce 
the almost sinusoidal profile of q^{kB) for the top band. 
The modes of the top band with small ks experience 
normal (as in free space) diffraction, while those with ks 
near tt have anomalous diffraction. This is in contrast 
to the system of coupled metallic slot waveguides, where 
the situation is the opposite due to e being negative in 
metals [9 . 

As we reduce 5, whilst keeping w = 220nm, we find 
that for TM modes, the gap below the top band shrinks 
and disappears at the Brilluoin zone edges. The gap clo- 
sure happens when s = {71-^63 + Cg — eskw)/eg/k = sq. 
This result is derived using the exact boundary condi- 
tions at the interfaces and is reproduced well while nu- 
merically solving Eq. (17) {sq ~ 125nm for w = 220nm). 
The critical value sq exists due to the Brewster angle 
condition, which gives zero reflection of the TM polar- 
ized wave from an interface. When the Brewster con- 
dition is satisfied, then the resonant Bragg scattering is 
canceled, the gap is closed and the periodic medium be- 
comes transparent [29l [30] . For s < sq the gap in q opens 
again, while the lower band continues to sink and soon 
disappears below the q^ = line. The band structures of 
TM waves and Bloch modes in this regime are shown in 
Fig-Htd). 

Crucially, at the instant when s < sq the geometry 
of the top band TM modes undergoes structural trans- 
formation. Specifically, the out-of-phase modes at the 
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FIG. 2: (Color online) Spectra of linear TE (a,c,e) and TM (b,d,f) waves in different periodic structures: (a),(b) w = 220nm, 
s — 500nm; (c),(d) w — 220nm, s — 50nm; (e), (f) w — 95nm, s — 50nm. Small plots show the mode profiles at the centre and 
edge of the of the Brillouin zone from the top band. The inset in (a) shows zoom of the top band. 



bottom of the band now cross zero not inside the low 
index material, but inside the high index one, cf. insets 
in Fig.|2|b) and (d). This transformation signals a qual- 
itative transition to the regime, where it is appropriate 
to consider our periodic medium as an array of coupled 
slot waveguides with sub- wavelength light localization in- 
side the low index slots. Indeed, the field intensity of the 
TM modes at the top of the band now peaks in-between 
the semiconductor layers and is strongly depressed inside 
them, see inset in Fig. [2|b). 

In contrast to the TM modes, the first band gap of 
the TE modes never shrinks to zero, and the structural 
transformation of the profiles of the top-band TE modes 
does not happen, cf. Figs. [2|a) and (c). The lower band 
of the TE spectrum sinks below the = line in a way 
similar to the TM case. 

Further reduction of either s or w (or both of them si- 
multaneously) leads to the edges of the transmission band 
sinking below the = cutoff for both the TE and TM 



modes, so that the out-of-phase modes corresponding to 
anomalous diffraction gradually disappear, see Figs. [2|e) 
and (f). The linear spectrum in this regime qualitatively 
reproduces that of a homogeneous medium, so that the 
description of the structure with an effective index ap- 
proach becomes relevant. The TE modes in this regime 
approach the ones of a homogeneous medium. E.g., the 
TE mode in the middle of the Brillouin zone tends to a 
constant, see inset in Fig. |2|e). At the same time the 
TM modes remain deeply modulated, with pronounced 
jumps at interfaces, see inset in Fig. [2|f). 

The excitation of a periodic medium with a narrow 
beam naturally results in a diffractive spreading. We can 
estimate the diffraction length Id of a beam with a radius 
d diS Id = k(P/\S\^ where the diffraction coefficient 



6 = 



dPq 
dk 



(19) 



Provided the initial excitation has a flat phase across the 
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beam, the 6 is calculated at /c^ = of the highest band. 
For s = 50nm and w = 220nm or w = 95nm the diffrac- 
tion length of the 500nm wide beam is approximately 
3/im and is roughly the same for TE and TM waves. 
However, for w = 500nm the for TM modes is still 
around 3/im, while for TE modes it is about an order of 
magnitude more. We note that the overlap of the TE 
modes through the evanescent fields inside the dielectric 
layers is smaller than the overlap of the TM modes. This 
is consistent with the difference of the diffraction lengths. 

Using the soliton concept we aim to demonstrate that 
preventing this fast spreading and achieving subwave- 
length localisation of light is possible. Below we are es- 
tablishing existence and study stability of solitons and, 
on this basis, are making some assumptions about their 
mobility. The propagation studies, which will allow full 
exploitation of soliton dynamics and functionality, re- 
quire development of a different set of numerical tools 
and are postponed for the future. 



IV. NUMERICAL METHODS FOR FINDING 
SOLITON SOLUTIONS AND DETERMINING 
THEIR STABILITY 



Soliton solutions of Eqs. ( 12 )-(|l4|) discussed in the next 



chapter have been found using the shooting method with 
zero boundary conditions at infinity. Since the nonlin- 
earities of silica and silicon are positive, the soliton prop- 
agation constant q has to be greater than the one for the 
linear waves, i.e. q > qnn. Here qnn is the propagation 
constant of the = linear mode of the top band. 

To characterize soliton solutions we use the power den- 
sity Pz defined as the x integrated z component of the 

time averaged Poynting vector ( S , 



f 

J — ( 



{S^) dx, S = £ xn 



(20) 



where (..) denotes time averaging. We plot as the 
function of the nonlinear phase shift induced by a soliton. 
The phase shift is defined as the difference between the 
total and the maximal linear wavenumbers 



(l) = k{q- qiin). 



(21) 



Substituting E{x^z) = {Eo{x) + e{x^ z)) ex.-p{iqkz)^ 
H{x, z) = {Ho{x) -\-h{x^ z)) ex.-p{iqkz) and linearizing the 
Maxwell equations for small e, h we find that the latter 
obey: 



idza 
kLb 



kAa + Kb . 
Ma. 



(22) 
(23) 



[^xi ^yi hy] , b — [(^zi 



q -l/(ceo) 

q l/(ceo) 

c€o{e + iyyy) q 
-ceo(e + J/xx) q 







" 


k = 











_ kceoUxz 







L = 


ceo(e 



ida: 
idx 



-l/(ceo) 




and 



27 

3 

27 



1^0 1^ + 2-Eio-E*o 



i = x,y,z 



(24) 
(25) 

(26) 

(27) 
(28) 



By noting that L can always be inverted and assuming 
a{x^z) = a(x) exp(— z/cAz), we arrive at the eigenvalue 
problem 



k Xd 



k^A^KL 



(29) 



Here Eq and Hq are the soliton solutions, a 



We approximate the derivatives in the matrices K and 
M by finite differences and are interested in the spatially 
localized a only. Then any A found with ImX > corre- 
sponds to an unstable perturbation exponentially grow- 
ing with the propagation coordinate z. 



V. NUMERICAL RESULTS 

For TM and TE waves we looked for and found two 
families of soliton solutions, which differ by the position 
of their center of symmetry. First, is the family of on-site 
solitons^ where the soliton center of symmetry is located 
in the middle of one of the semiconductor layers. So that, 
by saying site we assume a semiconductor layer. Second 
is the family of off-site solitons^ where the soliton center 
of symmetry is located between the semiconductor layers, 
i.e. within a slot filled with silica glass. 

On-site and off-site soliton profiles for w = 220nm and 
So < s = 500nm are shown in Fig. [3j One can see that 
for this relatively large separations, the TE and TM soli- 
tons have a similar field structure with the light mostly 
concentrated inside the semiconductor layers. Such field 
profiles can be well approximated by considering evanes- 
cently coupled semiconductor waveguides. The TM soli- 
ton is broader for the same power, since the modulus 
of the corresponding diffraction coefficient is by order of 
magnitude larger: Stm — —0.32 and Ste — —0.04. The 
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FIG. 3: (Color online) Intensity profiles of TE (black lines) 
and TM (red/grey lines) solitons: (a) on-site; (b) off-site. The 
power density for all solitons is fixed to = 10 GW/m, see 
the horizontal lines and markers A and B in Fig. [4] Grayscale 
background illustrates the underlying periodic structure: w = 
220nm, s 500nm. 
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FIG. 4: (Color online) Power density as a function of the 
phase shift for TE (a) and TM (b) solitons. Black (red/grey) 
lines correspond to the on-site (off-site) soliton. Solid 
(dashed) lines indicate stable (unstable) soliton branches. 
Structure parameters: w = 220nm, s = 500nm. Markers 
A and B correspond to the on-site and off-site solitons shown 
in Fig. [I 



on-site solitons are stable in this case, while the off-site 
ones are unstable. Note, that the on-site and off-site 
solitons in waveguide arrays described by the discrete 
nonlinear Schodinger equation have the same stability 
property [lOl [3j . Power as a function of is shown 
in Fig. [4j The on-site and off-site solitons make the same 
phase shift providing that the latter has a higher power. 
This is true for both the TE and TM families. 

As s approaches and becomes less than 5o, while w is 
kept fixed at 220nm, the TM solitons undergo a quali- 
tative change similar to the changes in the linear Bloch 
modes. The intensity peaks within the slots start to grow 
and prevail, so that the on-site TM solitons now have a 
two-peak structure, see Fig.[5|a,c), while the off-site ones 
have a single dominant peak, see Fig. [5jb,d). Note, that 
we are still in the regime, when the entire top band of 
Bloch modes has > 0. 

The power as the function of the phase shift (j) is 
shown in Figs. |6ja) and (b) for the TE and TM solitons, 
respectively. For both the TE and TM families, the on- 
site solitons still give a larger phase shift at a given power 
than the off-site solitons. The unstable eigenmodes of the 
off-site solitons from Fig.jsjd) are shown in Fig.[7|^a). In 
both the TE and TM cases, the instability is associated 
with the anti-symmetric eigenmode (known as "depin- 
ning mode" in tight-binding models ^Slj) and is therefore 
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FIG. 5: (Color online) Intensity profiles of TE (black lines) 
and TM (red/grey lines) solitons: (a,c) on-site; (b,d) off-site: 
w = 220nm, s = 50nm. The power density for the solitons 
shown is fixed to: = 33 GW/m (a,b) and P, = 110 GW/m 
(c,d), see the horizontal lines and markers A,B,C,D in Fig.[6] 
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FIG. 6: (Color online) Power density as a function of the 
phase shift for TE (a) and TM (b) solitons. Black (red/grey) 
lines correspond to the on-site (off-site) soliton. Solid 
(dashed) lines indicate stable (unstable) soliton branches. 
Structure parameters: w = 220nm, s = 50nm. Markers A-D 
correspond to the on-site and off-site solitons from Fig. [sja)- 
(d), respectively. 



expected to induce soliton motion across the structure, 
resulting in emission of dispersive waves and gradual con- 
vergence to a stable on-site soliton. 

We now reduce the width of the silicon layers to 
w = 95nm and enter the regime, when the linear diffrac- 
tion law starts to approach the limit of a homogeneous 
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FIG. 7: (Color online) Unstable eigenmodes for TE (black 
lines) and TM (red/grey lines) solitons: (a) off-site TE and 
TM solitons from Fig. [SJd); (b) off-site TE soliton from 
Fig.jsjd) and on-site TM soliton from Fig.jsjc). 
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FIG. 8: (Color online) Intensity profiles of TE (black lines) 
and TM (red/grey lines) solitons: (a,c) on-site; (b,d) off-site: 
w = 95nm, s = 50nm. Power density for the solitons shown 
is fixed to: = 50 GW/m (top row) and = 110 GW/m 
(bottom row), see the horizontal lines and markers A,B,C,D 
in Fig. [9] 
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FIG. 9: (Color online) Power density as a function of the 
phase shift for TE (a) and TM (b) solitons. Black (red/grey) 
lines correspond to the on-site (off-site) soliton. Solid 
(dashed) lines indicate stable (unstable) soliton branches. 
Structure parameters: w — 95nm, s — 50nm. Markers A-D 
correspond to the on-site and off-site solitons from Fig. [sja)- 
(d), respectively. 



material (i.e. when the edges of the top band sink below 
q"^ = line and anomalous diffraction regions disappear, 
see Figs.|2];e),(f)). In this regime, nothing new happens 
to the TE solitons and their shape feels only very lit- 
tle of the material inhomogeneity, see Fig. [8| At the 
same time, the structure of the TM solitons is still very 
imhomogeneous with sharp peaks within silica slots, see 
Fig.[8| As a significant fraction of light intensity in TM 
solitons is now concentrated inside slots with low n2 , TM 
solitons are noticeably broader than TE solitons at the 
same level of power, despite the corresponding diffrac- 
tion coefficient being larger by modulus for TE top band 
mode: 6tm — —0.3, Ste — —0.4. Powers as function of 
(p for this geometry are plotted in Fig. [9j While nothing 
qualitatively changes for TE solitons, on-site and off-site 
TM solitons exchange their roles: now they produce the 
same phase shift provided the former has higher power. 
Importantly, this power exchange is accompanied by the 
exchange in stability as well. The off-site soliton has 



become stable, while the on-site has got the instability 
with respect to the anti-symmetric linear eigenmode, see 
Fig. [rj^b). This result strongly suggests that the role of 
an elementary waveguide (or site) should be reconsidered 
for TM solitons in nanostructured periodic medium with 
the linear spectrum approaching the spectrum of the ho- 
mogeneous material. Indeed, since both linear and non- 
linear modes show strong light confinement inside slots, 
it is appropriate to consider a slot waveguide as an el- 
ementary structure. It means that if the coupled mode 
approach is to be developed for such geometries in the 
future, it should use the slot modes as the basis. 

We also note that stabilization of the off-site and desta- 
bilization of the on-site solitons with respect to the an- 
tisymmetric depinning perturbations happens in the dis- 
crete nonlinear Schrodinger equation if nonlinear cou- 
pling between the adjacent sites is included [32]. Remark- 
ably, in our case the coupling between the slot modes oc- 
curs through overlap of the evanescent fields inside the 
silicon, which is 100 times more nonlinear than the silica 
glass inside the slots. Similarly to other systems j3T| [32], 
the stability exchange between the on-site and off-site 
solitons is expected to be accompanied by an enhanced 
mobility of strongly inhomogeneous TM solitons. The 
possibility of such enhanced mobility is a problem opened 
for future investigation, which can reveal new possibili- 
ties for all-optical signal steering and manipulation in 
sub- wavelength structures. We have also checked for the 
existence of the cross-polarization instabilities, but have 
not found any for the examples presented above. This 
issue deserves further more detailed analysis. 



VI. SUMMARY 

Using the first principle nonlinear Maxwell equations 
we have developed and applied numerical techniques 
for finding solitons in periodic semiconductor-dielectric 
nanostructures and determining their stability. 

When the separation between the high index (semi- 
conductor) layers in a periodic structure is reduced be- 
yond the critical value (Brewster condition), the struc- 
ture of the antisymmetric TM modes at the edges of the 
highest Brilluoin zone undergoes qualitative transforma- 
tion. Namely the zeros of these modes now appear in 
the middle of the high index layers and not in the low 
index ones. This signals transition to the regime when 
a periodic nanostructure can be considered as an array 
of coupled slot waveguides for TM waves, rather than 
an array of coupled high index waveguides guiding by 
means of total internal reflection. Reducing the inter- 
face separation further one enters the regime where the 
diffraction law in a periodic nanostructure becomes sim- 
ilar to the diffraction in a homogeneous medium. The 
TE solitons in this quasihomogeneous limit feel little or 
no periodicity. However, the transverse profiles of the 
TM solitons reflect the subwavelength periodic structure 
and have dominant intensity peaks inside the low index 
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slots. In this regime the TM sohtons with the center of 
symmetry in the middle of a low index slot are stable 
while the ones centered on a high index layer are unsta- 
ble. This further facilitates the idea that the TM waves 



in periodic nanostructures can be considered as modes of 
the coupled slot waveguides and calls for developing of 
the corresponding analytic approaches. 
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